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At low values of external doping graphene displays a wealth of unconventional transport proper- 
ties. Perhaps most strikingly, it supports a robust 'metallic' regime, with universal conductance of 
the order of the conductance quantum. We here apply a combination of mean field and bosonization 
methods to explore the large scale transport properties of the system. We find that, irrespective of 
the doping level, disordered graphene is subject to common mechanisms of Anderson localization. 
However, at low doping a number of renormalization mechanisms conspire to protect the conduc- 
tivity of the system, to an extend that strong localization may not be seen even at temperatures 
much smaller than those underlying present experimental work. 
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The low energy physics of two-dimensional carbon, 
or graphene P, is governed by the presence of two 
generations of Dirac fermions. A conspiracy of these 
quasi-particles with various channels of impurity scat- 
tering leads to a wealth of intriguing transport phenom- 
ena which have attracted a lot of recent attention (see 
Refs. 00] for a list of early references.) Prominent exam- 
ples of the unorthodox conduction properties of graphene 
monolayers include a quantum Hall effect with uncon- 
ventional plateaux quantization Uxy = (2n+ l)e^//i, and 
the experimental observation of a minimal universal con- 
ductivity axx ^ /h at low concentration of dopands, 
iJ_F^T^^[l|,|2|. (T~^is the elastic scattering rate and 
Ep the Fermi energy, vanishing in the limit of zero exter- 
nal doping.) While many of graphene's transport prop- 
erties have been successfully addressed theoretically, a 
number of important questions remain open. Specifically, 
the nature of the low energy regime Ep ^ — where 
the deviations from the behavior of conventional metals 
are most pronounced — is not fully understood. Relying 
on calculations to lowest order in perturbation theory in 
the disorder concentration the experimental results 
have been interpreted in terms of some delocalized, 
genuinely metallic phase at the band center. However, 
in the very low energy regime, where the conductivity 
begins to level of to a disorder independent value, per- 
turbation theory breaks down and does not really help 
to elucidate the physics of the system. 

Building on a combination of mean field analysis and 
bosonized perturbation around the clean Dirac fermion 
fixed point, Ep = 0, we here propose an effective low 
energy field theory of system. Refiecting the large scale 
symmetries of the problem |^, the relevant field theory 
at strong doping {EpT 3> 1) is a standard nonlinear a- 
model of either orthogonal or unitary symmetry (depend- 
ing on whether time reversal invariance is broken or not.) 
The bare coupling constant of the theory (cf. Eq. ^ be- 
low) Gxx EpT ^ 1 is large, implying that localization 
corrections to the Drude conductivity are comparatively 
weak. Our main finding is that the weak doping regime, 



EpT <C 1, too, is described by the nonlinear cr-model. 
This result implies that at all doping levels disordered 
graphene is subject to common mechanisms of Ander- 
son localization. Specifically, the smallness of the bare 
coupling constant controlling the theory at EpT < 1 
— loosely to be identified with the disorder independent 
'minimal universal conductivity' a^x = C(l) of the sys- 
tem — signals that localization and fluctuation phenom- 
ena will be especially pronounced at small doping. Below 
we will discuss a number of mechanisms which may ex- 
plain why these fluctuations have not yet been observed 
experimentally. We conclude by briefly considering the 
Quantum Hall physics of the system. 

In the absence of disorder, graphene's two generations 
of low energy Dirac fermions are described by the Hamil- 
tonian Q H^s = —ivdxOx^T'i—ivdxCry, where v ^ ta, a is 
the lattice spacing, the Pauli matrices Ui act in the space 
of the two sublattices A and B of the system, and the 
Pauli matrices r, act in two component nodal space. The 
off-diagonality of the Hamiltonian in sublattice space im- 
plies chirality, [H,az\+ = 0. 

In the presence of all possible channels of disorder scat- 
tering, H generalizes to (outer block structure in nodal 
space, 2x2 sub-blocks in sublattice space), 
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where d — dx + idy, and Vc and Wc, respectively, de- 
scribe the soft (intra-node) and the hard (inter-node) 
scattering due to defect sites, next nearest neighbor 
hopping, or disorder in adjacent substrates. Diagonal 
in sublattice space, these scattering amplitudes violate 
the chiral symmetry of the clean system. Similarly, 
the chiral symmetry preserving matrix elements V and 
, W represent the intra- and inter-node ampli- 
tudes, respectively, of defective bond hopping. 

Referring for a much more comprehensve discussion, 
to Ref. , we briefly comment on the symmetries of the 



2 



problem: (i) if all hopping matrix elements of the original 
problem are real (complex), the low energy Hamiltonian 
does (not) obey the symmetry relation H = tiH^^ti, thus 
falling into the Wigner-Dyson orthogonal (unitary) sym- 
metry class. However, (ii) depending on the microscopic 
realization of the disorder, some of the amplitudes en- 
tering the bare Hamiltonian may be much smaller (even 
negligibly small at the length scales probed in current ex- 
periment) than others. If this happens the approximate 
realization of symmetries beyond Wigner-Dyson orthog- 
onal/unitary may result in various transport anomalies, 
including weak antilocalization phenomena at intermedi- 
ate length scales However, not protected by 
rigid exclusion principles, these structures are transient 
in nature and will not affect the physics at the largest 
distance scales. As we are primarily interested in the lat- 
ter, we will assume the presence of all disorder channels 
throughout. 

Transport properties of the system are described 
by disorder averaged products of retarded and ad- 
vanced Green functions {{Ep — H + iO)xi,x2(^P — H — 
*0)x3^,x4 • ■ ■ )dis- To compute such correlation functions 
we consider the replica coherent state field integral (gen- 
eralization to a supersymmetry formulation is straight- 
forward) 
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where ip — {ips,a,!y{^)} is a (2 x i? x 4)-component 
Grassmann field, s — 1,2 is a two component index 
distinguishing between retarded and advanced indices, 
a = 1, . . . , i? is a replica index, v — 1, . . . , 4 enumerates 
sublattice and nodal sectors and trf' = (Jss'(— )'*^^ is a 
Pauli matrix causing infinitesimal symmetry breaking in 
advanced/retarded space. 

Beginning with the strong doping regime^ EpT ^ 1, we 
next outline the construction of the low energy field the- 
ory describing the large scale behavior of the system (i.e., 
at length scales larger than the elastic mean free path, 
I.) The construction of the theory at large doping follows 
standard procedures (cf. Refs. [ulllJI for review) and we 
restrict ourselves to a brief recapitulation of the central 
construction steps: the action in (j^J is invariant under 
uniform unitary rotations "0 ~^ Uip, where U £ U(2i?). 
Averaging over all scattering channels in |^ generates 
a mean field fermion self energy i^crf jScrf ', break- 
ing the symmetry down to U(i?) x U(i?). Along with 
this spontaneous symmetry breaking, Goldstone modes 
Q e V{2R)/V{R) X U(i?) appear. The effective ac- 
tion S[Q] controlling the low energy partition function 
Zcs = J DQ exp(— 5[(3]) of these modes is fixed by not- 
ing that the field space admits a (unique) parity invariant 
second order derivative operator, — J (Pxti (d^Qdf^Q). 
Gauge arguments (but also the microscopic construction) 
state that the bare coupling of this operator is (1/8) x the 



dimensionless Drude longitudinal conductivity, ctxx, i-e. 

S[Q] = J d^xtv {d^Qd,,Q), (3) 

which is the standard nonlinear a-model action of time 
reversal invariance broken weakly disordered metallic 
systems. (In the time reversal invariant case, Q £ 
Sp(4i?)/Sp(2i?) X Sp(2i?) and a^^/S cr^^/lQ.) Specif- 
ically, graphene's Drude conductivity is given by 
a XX ^'^Z^^ 2tEp. Referring for an in depth discussion of 
the localization properties of the theory © to Ref. [l^, 
we note that at intermediate length scales L > I this 
bare value will be modified by weak localization correc- 
tions (for a comprehensive discussion of weak localization 
effects in graphene at Eft > 1 and different channels 
of disorder scattering, see Ref. 0), before strong An- 
derson localization sets in at exponentially large scales 
L ^ = Zexp(const. x cr^^) (unitary symmetry) or 
^ = Zexp (const, x axx) (orthogonal symmetry.) How- 
ever, before discussing the phenomenology of the system 
any further, let us explore what happens in the strong 
disorder, or weak doping regime, EpT < 1. 

While the mean-field type derivation of the action Q 
essentially relied on the presumed largeness of the pa- 
rameter EpT, changes in the bare disorder concentration 
are not expected to change the basic symmetries of the 
system. One may, thus, suspect that the low energy the- 
ory retains its integrity as we push on into the strong 
disorder regime EpT < 1. In the following we will apply 
bosonization techniques (similar to those previously used 
in the context of the d-wave superconductor [ij, [ij] ) to 
confirm that this conjecture is correct. Readers not in- 
terested in the technicalities of this construction may di- 
rectly advance to the discussion two paragraphs further 
down. 

Our starting point is Wittens seminal result ac- 
cording to which a system of N free Dirac fermions may 
be equivalently represented in terms of a level 1 Wess- 
Zumino-Witten (WZW) model with target manifold 
U{N): J d^x + ii-d'^-) ^ S{M\, where V'± are 

A^-component Grassmann fields and M S U(2i?) is an 
N X N unitary matrix field with action S[M] = Sq[M] + 
j^T[M]. Here, So[M] = ^/ d^x tr (9^ MS^M"!) and 
r[M] is the topological WZW-action whose detailed 
structure will not be of relevance throughout 0| . Specif- 
ically, the replica theory of clean graphene with its two 
generations {K and K') of 2R Dirac fermions is repre- 
sented in terms of two matrix fields M,M' G U(2i?), 
with decoupled action S[M,M'] = S[M] + S'[M']. Here 
S[M] = So[M] + ^V[M] and S'[M'] = So[M'] - 
-p^T[M'] where the relative sign in front of the WZW- 
actions signals the different orientation (or parity) of the 
two low energy islands. The role of disorder scattering 
may now be explored by averaging the exponentiated im- 
purity matrix elements in the prototypical action |(2Jl over 
the disorder distribution functions (here assumed to be 
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Gaussian correlated.) This results in the appearance of procedure on the inter-node chiral symmetry preserving 
four-fermion operators, which may be subsequently be scattering amplitudes, we find 
represented in bosonic language [l5|- Exemplifying this 

I 

\ I dis. 

I 



where the constant "f^^ is a measure of strength of 
the disorder, a is the small-distance cutoff of the the- 
ory (the lattice spacing, say), and the correspondence 
ip-aip+b CL~^Mab was uscd. In a similar man- 
ner, the scattering channels Va.b generates operators 
-7Aa-2 J [tr (M + M'-^)] '\ 

To understand the consequences of the disorder gen- 
erated modification of the theory, we note that the 
two operators above are examples of marginally rele- 
vant perturbations at the clean Dirac fermion fixed point. 
For initially weak disorder, the con- 
stants J A and grow to values of 0(1) only at ex- 
ponentially large length scales / ^ aexp(const. x w^/70), 
where 70 represents the bare values of the coupling. The 
crossover scale, I, defines the elastic mean free path of 
the theory. At larger scales, the inter-node scattering 
operator tr (MM') enforces locking M = M' of the for- 
merly independent nodes while node and sublattice di- 
agonal scattering, [tr(M -I- M')]^^, implies a field reduc- 
tion M ^ Q G U(2i?)/U(i?) X \J{R) down to the mani- 
fold discussed earlier in the context of the weak coupling 
regime. Substituting the reduced fields into the action 
S[M,M'] — > S[Q,Q] (For the legitimacy of naive mass- 
less field substitution within a strong coupling context, 
see Ref. |i3) we obtain Eq. © at axx — 4/7r. (For time 
reversal invariant systems, a similar construction leads to 
the orthogonal variant of the theory.) 

We conclude that the theory describing electronic 
transport in graphene at low doping, EpT <^ 1, is a non- 
linear CT-model with disorder independent bare coupling 
constant i/ir. However, unlike with the weak coupling 
situation discussed above, (a), this description becomes 
vaHd only at length scales I ~ aexp(const. x Vp/^o), 
exponentially large in the bare inverse disorder concen- 
tration. At smaller length scales, the dynamics is es- 
sentially ballistic. (Within the framework of a large N 
mean field approach, this strong disorder sensitivity of 
the Dirac fermion scattering mean free path, and the rel- 
evancy of the cr-model at larger scales was first observed 
in Ref. 0.) In contrast, the mean free path I 
at strong doping EpT ^ 1 depends much weaker on 
the disorder concentration and the crossover to a mul- 
tiple scattering regime occurs at smaller length scales. 
Second, (b) the smallness of the coupling constant A/tt 
implies that the bare value of the conductivity will be 



subject to strong renormalization. More precisely, large 
fluctuations of the Q-fields (describing the proliferation 
of the particle-holc modes causing Anderson localization) 
will lead to an exponential decay of correlation functions 
at length scales ^ = J\fl, where A/" is a numerical fac- 
tor. Accordingly, the conductivity of an infinite graphene 
sheet axx{T = 0) = 0, while the conductance of a sam- 
ple with extension L ^ I will be exponentially small in 
exp(— const, x L/^). However, as first pointed out in 
the value min (L, ^0) beyond which the conductivity is 
strongly diminished may be orders of magnitude larger 
than /: naive application of a 1-loop RG to the strong 
coupling problem ^ suggests that for systems with bro- 
ken time reversal invariance, localization sets in at scales 
^ ^ e^^l. In time reversal invariant systems ^ ~ e*Z, still 
about two orders of magnitude larger than /. 

How can these results be reconciled with the experi- 
mental observation of a universal value axx,cxp — 4, sig- 
nificantly larger than even the bare (Drude) conductivity 
of the clean Dirac system? At this stage, the numerical 
discrepancy to the Drude value is not fully understood 
(see, however, Ref. ^3-) general robustness 

of the minimal conductivity, (i) elastic mean free path 
renormalization at small doping and (ii) the largeness of 
the sigma model correlation length even at strong cou- 
pling conspire to render the length scale beyond which 
strong localization is seen large. Finally, (iii), we repeat 
that the above discussion applied to a regime where all 
disorder scattering channels are of comparable strength. 
On the other hand, the experimentally observed linear 
dependence of the conductivity on the carrier density 
suggests ^3 dominance of the soft (node diagonal) po- 
tential scattering, Va = Vb channel. If only this chan- 
nel were operational, the relevant field theory would be 
two independent copies (the nodes) of sigma models of 
symplectic symmetry (technically, Witten's theory sub- 
ject to the constrait M = Q, M' — Q'), each gener- 
ating weak aniilocalization. However, beyond a certain 
crossover scale, internode scattering will become effec- 
tual, the field theory gets locked to a single copy of a 
theory of orthogonal symmetry and standard Anderson 
localization emerges. We thus conclude that graphene is 
subject to ordinary mechanisms of Anderson localization. 
Although, a conspiracy of points (a-d) above may imply 
that the observed metallic regime is extraordinarily ro- 
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bust, weak localization precursors of strong localization 
will likely be observed at temperatures lower than the 
present 10 K. 

Slightly modifying a line of arguments originally due 
to Pruisken|19l|. we finally touch upon the quantum 
Hall regime. On very general grounds, i.e. regard- 
less of whether a magnetic fields is present or not, 
the theory contains a first order derivative operator, 
5(1) [T] = J d^x {j^i^)piEF))tT{al'T-^d^T), where 
jp(x) is the local current density operator, p{Ep) the 
(disorder averaged) density of states at the Fermi energy 
and {■■■)q the trace over single particle Hilbert space. 
While the bulk expectation value {jfj,p)q generally van- 
ishes, the presence of quantum Hall edges couples 5'^^^ 
to the theory. Substitution of the edge state system of 
graphene j20lEH into S^^''[T] obtains the boundary oper- 
ator ^(i) [T] = -v§ds (tr(crf T-i^^T), where s parame- 
terizes the system boundary and = 2n + 1 is the num- 
ber of Landau levels below the Fermi energy. Applying 
Stokes theorem to convert S'^'^^ to a bulk integral, adding 
the background action (Q, and using that v = Oxyl'^ for 
Fermi energies intermediate between Landau levels, we 
obtain 

S\Q\ = J d'^xtr [a^^df^Qdi^Q - crxy(^iiuQd^Qdf,Q] . 

(4) 

For high Landau levels, Q is but a variant of Pruiskens 
quantum Hall action; the bare longitudinal conductivity 
axx is large and the structure of graphene's edge spec- 
trum merely enters through the odd-integer quantization 
of the Hall conductivity. With the lowest Landau level 
(LLL), V = 1, the situation is more interesting. Refer- 
ring for the technicalities of bosonization in the presence 
of magnetic fields and system boundaries to Refs. 
we note that the applicability of the action |0J to the LLL 
must be taken with a grain of salt: for one thing, the 
bosonization approach perturbs around the clean, non- 
magnetic limit of the theory. Quantities such as the dis- 
order averaged broadened density of states of the LLL 
are likely inaccessible in terms of (renormalized) pertur- 
bation theory around that point. Second, and unlike with 
the situation in ordinary metals, the conductivity tensor 
of the LLL {(Jxx,cFxy) = 2 x (2/7r, 1) is numerically close 
to the prospected value of the quantum Hall transition 
point (cr*^,, (J%y)- This is important inasmuch as the QHE 
transition is arguably |23ll24| not described by Pruiskens 
theory. In practice, this may imply that strong field fluc- 
tuations (to be expected on account of the smallness of 
the coupling constants) readily compromise the integrity 
of the LLL action Q). At any rate, the identity of the 
proper theory of the quantum Hall transition remains as 
unknown as it is in normal metals. 

Summarizing, we have constructed the low energy the- 
ory of disordered graphene at large length scales. In its 
low doping regime, EpT < 1, the system is described 
by a nonlinear cr~model at strong coupling. Derived by 



bosonization methods, this model consistently matches 
the weak coupling mean field approach, valid at large 
doping EpT ^ 1. Our main finding is that the 'uni- 
versal minimal conductivity' in graphene will be subject 
to conventional mechanisms of quantum interference at 
temperatures lower than those underlying current exper- 
iment. However, depending on the symmetries actually 
realized in the system, the temperatures below which the 
conductance is strongly diminished by localization may 
be extremely low. 
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I thank E. Fradkin for pointing out to me the pro- 
found applied relevance mean free path renormaliza- 
tion likely has to the physics of the low doping regime. 
This work was supported by Transregio SFB 12 of 
the Deutsche Forschungsgemeinschaft. Note added: af- 
ter completion of this work, I became aware of the 
closely related manuscript I.L. Aleiner and K.B. Efetov, 
.cond-mat/0607200| 



[1 
[2: 

[s: 
[4: 

[6; 

[7: 

[9 

[10: 
[11 

[12 

[13 

[14 
[15 
[16 



245411 



[17] 
[18] 
[19] 
[20] 

[21] 
[22] 

[23] 
[24] 



K. Novoselov et al, Nature 438, 197 (2005). 
Y. Zhang et al, Science 438, 201 (2005). 
V. Gusynin and S. Sharapov, Phys. Rev. B 95, 
(2005). 

T. Ando, J. Phys. Soc. Jpn. 74, 777 (2005). 

V. Gusynin and S. Shaparov, Phys. Rev. Lett. 95, 146801 

(2005). 

H. Suzum-a and T. Ando, Phys. Rev. Lett. 89, 266603 
(2002). 

D. Khveshchenko, cond-mat/0604180. 
McCann et al, cond-mat/0604015. 
Fradkin, Phys. Rev. B 33, 3263 (1986). 
Nomura and A. MacDonald, cond-mat/0606589. 
Altland and B. Simons, Condensed Matter Field The- 
ory (Cambridge University Press, 2006). 

K. Efetov, Supersymmetry in Disorder and Chaos (Cam- 
bridge University Press, 1997). 

A. Nersesyan et al, Phys. Rev. Lett. 77, 2628 (1994); 

Nucl. Phys. B 438, 561 (1995). 

A. Altland et al., Phys. Rep. 359, 283 (2002). 

E. Witten, Comm. Math. Phys. 92, 455 (1984). 
Compactifying the two-dimensional graphene sheet at in- 
finity to a two-sphere, S, the WZW-action is given by 
r[M] = iJgd^xe^^^tT{M-^d^MM-''d^MM-^daM), 
where B is the interior of S and M — M(x,s) is a ho- 
motopic extension of M{x) from S to B, subject to the 
condition M(x, 1) = M[x) and M[x,Q) = 1. 
A.W.W. Ludwig et al, Phys.Rev.B 50, 7526 (1994). 

L Affleck and F. Haldane, Phys. Rev. B 36, 5291 (1987). 
A. Pruisken, Nucl. Phys. B 235, 277 (1984). 
A. H. Castro Neto and N. P. F. Guinea, Phys. Rev. B 
73, 205408 (2006). 

L. Brey and H. Fertig, Phys. Rev. B 73, 195408 (2006). 
D. Delphenich and J. Schechter, Int. J. Mod. Phys. A 12, 
5305 (1997). 

M. Zirnbauer, hep-th/9905054. 

M. Bhaseen, Nucl. Phys. B 580, 688 (2000). 



